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Abstract 

o 

We propose a generalized feedback model and compressive sensing based opportunistic feedback 
schemes for feedback resource reduction in MIMO Broadcast Channels under the assumption that both 
uplink and downlink channels undergo block Rayleigh fading. Feedback resources are shared and are 
opportunistically accessed by users who are strong, i.e. users whose channel quality information is above 
a certain fixed threshold. Strong users send same feedback information on all shared channels. They are 
identified by the base station via compressive sensing. Both analog and digital feedbacks are considered. 
The proposed analog & digital opportunistic feedback schemes are shown to achieve the same sum-rate 
throughput as that achieved by dedicated feedback schemes, but with feedback channels growing only 
logarithmically with number of users. Moreover, there is also a reduction in the feedback load. In the 
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analog feedback case, we show that the propose scheme reduces the feedback noise which eventually 
results in better throughput, whereas in the digital feedback case the proposed scheme in a noisy scenario 
achieves almost the throughput obtained in a noiseless dedicated feedback scenario. We also show that 
for a fixed given budget of feedback bits, there exist a trade-off between the number of shared channels 
and thresholds accuracy of the feedback SINR. 



I. Introduction 

Recently, it has been shown that dirty paper coding (DPC) achieves the sum-rate throughput of the 
multiple-input multiple-output (MIMO) broadcast channel [1], [2]. However, it requires a great deal of 
feedback as the transmitter needs perfect channel state information for all users and is computationally 
expensive [3]. Since then, many works have attempted to achieve the same sum-rate throughput with 
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imperfect channel state information (reduced feedback load). This was done by applying opportunistic 
communication in the forward link [4]- [7]. 

By reviewing the feedback protocols suggested in literature, one can note that generally the following 
three components are fed back and user selection is based on either one or a combination of these 
components [4]- [9]: 

1) Channel Direction Information (CDI), e.g., beam index(BI), quantized channel index (QCI)... etc. 

2) Channel Quality Information (CQI), e.g., SNR, SINR, channel norm etc. 

3) User identity (ID) 

Feedback schemes can be differentiated according to whether the feedback is analog or digital. It is 
termed digital if the feedback involves only digital data (integer or bits) and analog otherwise. Feedback 
schemes can also be classified as either opportunistic or non-opportunistic. 

While the forward link is opportunistic in nature, almost all feedback schemes are non-opportunistic. 
Here, each user has a dedicated feedback channel. For example, in random beamforming (RBF) scheme 
proposed by Sharif and Hassibi [4], the users are differentiated according to the channel direction (i.e. 
what beam direction is the channel mostly aligned with) and accordingly users feedback to Base Station 
(BS) the SINR corresponding to that direction only. So, each user feedback one integer and one real 
number. In order to reduce the feedback load further Diaz et. al. in [7] propose a threshold based RBF. 
Here, instead of feeding back the SINR for the best beam for each receive antenna (one real plus one 
integer numbers), the user only transmits one bit to the BS, indicating whether or not the SINR on a pre- 
selected beam for any receive antenna is above a given threshold. The scheme is repeated for each beam. 
Since interaction or cooperation among the competing users is not allowed, hence defying opportunism, 
there is a linear increase in the feedback resources (or channels) with the number of users [5], [6]. Even 
if thresholding is applied, there is no reduction in the number of feedback channels. This is because the 
channels are reserved even when users are not sending any feedback information. 

Recently, some works have started to consider opportunistic feedback schemes where feedback re- 
sources are shared and are opportunistically accessed by strong users i.e. users whose CQI is above the 
given thresholds. Thus, in [8], Tang et. al. propose a feedback scheme with fixed number of feedback 
slots (channels) that are randomly accessed by strong users. In every slot, each strong user independently 
attempts to send back to the BS a data package containing its user identity (ID) with a probability. If two 
or more users feedback in the same slot, collision occurs and the feedback in that slot is discarded. In the 
case when multiple users successfully feeding back, the BS randomly selects one of the successful users. 
Although the scheme requires only an integer feedback per slot, it is suboptimal as the user is selected 
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randomly. The scheme was only proposed for the single-input single-output (SISO) case. In [9], Rajiv et. 
al. propose a feedback scheme based on random access slots for the MIMO case that requires only user 
identity feedback (also an integer feedback per slot). In this scheme, time is divided into slots that are 
equally divided among the beams (CDI). Each slot then corresponds to a pre-determined threshold. Thus, 
if a user's CQI (e.g. SINR) on a particular beam exceeds the threshold corresponding to a particular slot, 
that user would feedback on that slot. 

Both of the two schemes above require accurate timing-synchronization to avoid collisions, which is 
difficult to achieve in practice. Moreover, feedback to the BS is successful if there are no collisions, 
i.e., only one user is attempting to feed back in a slot. In addition, the two schemes only work for 
digital feedback but not when the designer is interested in analog feedback. In all the feedback schemes 
discussed above, the feedback links were assumed ideal when the forward links were subject to both 
fading and noise. This asymmetry in the way the two links are treated is unrealistic. 

In this paper, we consider a broadcast scenario where the forward and the feedback links are symmetric 
in that they are both i) non-ideal and ii) opportunistic or shared. Thus, both links undergo Rayleigh 
fading and are subject to additive Gaussian noise. Moreover, the channels in both links are shared and 
are opportunistic in the sense feedback channels are dominated by strong users. Finally, the feedback 
links can be used for both analog and digital feedback. 

The paper proposes a generalized feedback model and compressive sensing (CS) [10]- [14] based op- 
portunistic feedback protocols for feedback resource reduction. Just as in all existing feedback techniques, 
a number of channel directions or beam is first determined. For each direction, the number of feedback 
channels is fixed and strong users fedback their CQI information on all feedback channels. In the analog 
feedback case, each strong user feds back CQI value whereas in the digital feedback case, each strong 
user feds back "1" if his CQI is above a particular threshold and remains silent otherwise. This creates an 
undetermined system of equations in a sparse vector of users. We use the emerging compressive sensing 
technique to identify users who have fed back and to estimate the fedback CQI. Users with higher 
value CQI have a stronger chance of being recovered. The results obtained via compressive sensing are 
refined using least-squares. As the feedback links are noisy, so the BS backs off on the noisy CQI based 
on the variance of the noise. We obtain the optimum back off on the noisy CQI that maximizes the 
throughput. A user among strong users is selected (strongest in the analog feedback case & randomly in 
the digital feedback case). The scheme is repeated for each channel direction. Although we have used 
SINR feedback, the proposed schemes can work with any kind of CQI (e.g. SNR). It is important to note 
that our scheme is less sensitive to timing-synchronization errors, as the scheme will be affected only 
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if out of synchronization user is selected for a particular channel direction (the probability of which is 
low). 

The remainder of the paper is organized as follows. In section [TTJ generalized feedback model is 
introduced. In section [In] we discuss the proposed feedback strategy. In section JVl we present the sum- 
rate throughput obtained by the proposed schemes in the RBF case. In section [V] performance evaluation 
of the proposed feedback schemes is presented. Feedback channel training is discussed in section [VT] 
followed by numerical results and conclusions in sections IVIII and IVIIII respectively. 

Notation: We use bold upper and lower case letters for matrices and vectors, respectively. A T , A* 
and At refers to Transpose, Hermitian conjugate and pseudo-inverse of A respectively. E x [-] denotes 
the expectation operator w.r.t. x, and P[ ] is the probability of the given event. The natural logarithm is 
referred to as log(-), while the base 2 logarithm is denoted as log 2 (-). f(x) = 0(g(x)) is equivalent to 
f(x) = cg(x) where c is a constant. \A\ denotes the size of a set A. 



A. Downlink Transmission Model 

We consider a single cell multi-antenna broadcast channel with p antennas at the base station (transmit- 
ter) and n users (receivers) each having one antenna. The channel is described by a propagation matrix 
which is constant during the coherence interval and is known completely at the receiver. Let u £ C pxl 
be the transmit symbol vector and let X{ be the received signal by the i-th user, the received signal by 
the i-th user can then be written as 



where hj £ C p is the channel gain vector between the transmitter and the user, and Wi is the additive 
noise. The entries of hj and Wi are i.i.d. complex Gaussian with zero mean and unit variance, CM (0, 1). 
Moreover, u satisfies an average transmit power constraint E{u*u} = 1 and pi is the SNR of the i-th 
user. A homogeneous network is considered, in which all users have the same SNR, i.e. pi = p = P/p 
for i = 1, . . . , n, where P is the total power available at the transmitter assuming that the noise power 
is unity. We also assume that the number of mobiles is greater than or equal to the number of transmit 
antennas, i.e., n > p, and that the BS selects p out of n users to transmit to. 

B. Generic Multi-antenna Feedback Channel 

We present here a general model for the multiuser feedback channel with r feedback channels (possibly 
shared) among n users, in which users report channel quality information (CQI) to the base station in 



II. System Model 




i = 1, . . . , n 



(1) 
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order to exploit multiuser diversity. As we shall soon see, this model encompasses the existing feedback 
models. The feedback channels are described by a propagation matrix A which is constant during the 
coherence interval and is assumed to be perfectly known at the BS (receiver), and are to be independent 
of the downlink channel. Let v S C raxl be transmit feedback vector and let t/i be the signal received 
via the i-th feedback channel. The signal received through the i-th feedback channel is mathematically 
described as 
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or equivalently 

y = Av + w (2) 

where r < n and ay represents the (generally complex valued) gain of the i-th channel for the j-th user. 
Note that in contrast to the majority of existing feedback reduction techniques, a noisy feedback channel 
is assumed. The entries of w represent the additive noise and are assumed to be i.i.d. complex Gaussian 
with zero mean and variance a 2 , CM (0, a 2 ). 

If no fading is considered (A is deterministic), all entries of A are equal to a constant, whereas if the 
feedback channels undergo block Rayleigh fading, A remains constant during the coherence interval and 
its entries are i.i.d. complex Gaussian with zero mean and unit variance, ~ CAf (0, 1), assumed to be 
known perfectly at the BS via feedback channel training (discussed in Section IVD). 

We summarize in Table H how our feedback model © applies to the feedback models (opportunistic or 
not) suggested in literature [4]- [9]. Thus, in the non-opportunistic feedback model, each user is allocated 
its own feedback channel and the uplink channel matrix A becomes diagonal and of size n (equal to 
the number of users). For the opportunistic models proposed in [8] by Tang et. al., the feedback channel 
matrix A becomes diagonal of size r x r, where r is the number of feedback slots and is less than 
n. v represents feedback data in each slot, and when a collision in a particular slot takes place, the 
corresponding entry of v is not valid. The same model holds for [9] except that in this scheme r is not 
fixed but varies randomly. Also, r may not necessarily be less than n. In all these schemes, the additive 
noise w is set to zero. 

In this paper, we take a more general approach and consider a contention-based feedback protocol, 
which assigns independent multi-access contention channels for CQI reporting in which the feedback 
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process should itself be a filter that selects strongest users. There are different ways to interpret the 
system of equations (cf. (0). One possibility is to assume that each user is equipped with one antenna 
and the BS is equipped with r antennas. In this case aij represents the gain from the j-th user to the 
i-th antenna (spatial feedback channels). Another possibility is to assume that each single-antenna user 
is going to feedback the same information over r frequency bands shared with the other users. Thus, 
represents the gain of the j-th user in the i-th band (frequency feedback channels). 



Before we discuss the proposed feedback strategy, we present important compressive sensing results 
used in our work. A short introduction to compressive sensing is given in the Appendix. 

A. Sparsity Pattern Recovery Results 

Compressive sensing refers to the recovery of the sparsity pattern S (with |5| = s) of signal v G 1" 
accurately from limited measurements 



Two approaches for recovering the sparsity pattern in the noisy setting (cf. (0) are discussed here, the 
only exception being that these results are derived for the case when the entries of A and w are i.i.d. 



1) Sparsity pattern recovery using LASSO [13]: A recent paper by Wainwright [13] shows that it is 
surprisingly possible to recover the sparsity pattern of signals accurately from limited measurements in 
a noisy setting using LASSO which is l\ -constrained quadratic program (QP). The LASSO gives the 
estimated v. 



where a > is a user-defined regularization parameter. 

The number of measurements (or channels) required for successful sparsity pattern recovery using 
LASSO must satisfy r\ = cislog(n — s), where c\ is a constant. 

2) Sparsity pattern recovery using Maximum Correlation [14]: Similar results for sparsity pattern 
recovery from limited measurements in a noisy setting using maximum correlation, a much simpler 
method compared to LASSO, are derived by Fletcher et. al. in [14]. Maximum Correlation estimate is 
defined as the indices corresponding to the s largest values of /, where / is defined as follows 



III. Proposed Feedback Strategy 



S = {ie{l,...,n}\vi^0} 



(3) 



real Gaussians i.e., ~ N (0, 1) and Wi ~ 




(4) 
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The number of measurements (or channels) required for successful sparsity pattern recovery using 
maximum correlation must satisfy T2 = C2slog(n — s), where C2 is a constant. 

3 ) Refining CS results using Least Squares: In our framework, we propose to refine the results obtained 
via compressive sensing through least squares (LS) using the following procedure: once the sparsity 
pattern S is known, we can form matrix A5 with columns of A corresponding to S and estimate v as 
vls = A^y 

B. General Strategy of Using Compressive Sensing for Feedback 

Any feedback scheme has two components, a direction component and a magnitude component. The 
transmitter usually has certain pre-determined directions for which it seeks user feedback. Thus, the BS 
announces that it is seeking feedback for a particular direction. At this instant, the users whose channels 
lie at or are close to this direction, feedback their CQI (SNR, SINR, channel strength etc.). Now a limited 
number of users will feedback on the set of shared feedback channels according to input/output equation 



Thus the vector v in (J2J) is sparse with sparsity level determined by the number of users who feedback. 
CS can now be used to recover the sparsity pattern of v [13]- [14] (i.e. which user prefer that particular 
direction) and could also recover the vector v itself [13]) (i.e. users' feedback CQI). Moreover, the larger 
the value of particular CQI, the higher the chances of its recovery. Another factor that enhances the level 
of recovery is how sparse the vector v as compared to the number of feedback channels available. We 
need at least one strong user (i.e. s > 1) for each beam or direction in order to achieve full multiplexing 
gain which implies that small values of s are sufficient. To reduce the number of users who feedback s, 
we pursue a thresholding strategy where the user will feedback if his CQI is greater than a threshold £ 
to be determined. 

Now consider a particular beam (CDI) (all beams will behave in an identical manner as the users are 
i.i.d. and the beams are equi-powered). Noting that the users' CQI are i.i.d., we can choose ( to produce 
a sparsity level s. This happen by requiring that 



where F(Q or u is the complementary cumulative distribution function (CCDF) of CQI (SINR) defined 



F(C) = are max ( 
«e(o,i) \s 




(5) 



as: F(C) = P[SINR > C] 



cx p(-C/p) 
(i+C)"- 1 



C>o. 




Proof: Let ip = (™)u s (l — u) n s . Differentiating tp w.r.t u and setting the derivative to 0, and solving 
for u yields u = s/n. Thus, F((^) = s/n, or £ = F~ 1 (s/n). 
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C. Feedback Protocol for the Analog Feedback Case 

In the analog feedback scenario, users above threshold feed back their analog CQI value. The CS 
strategy then allows the BS to recover all users who transmitted their CQI. This off course will be true 
provided that the number of users who feedback is less than or equal to s. Here, we assume that the 
probability the a user is strongest for more than one beam is negligible as the number of users are 
relatively much larger than the number of beams. It has been shown in [4] that this is a valid assumption 
under these conditions. The steps of the proposed compressive sensing based opportunistic feedback 
protocol are as follows: 

1) Threshold Determination: BS decides on thresholding level ( based on the sparsity level that can 
be recovered. 

2) User Feedback: Repeat the following steps for each beam. 

• CQI Determination: Each user determines his best beam (corresponding to the highest CQI 
value). 

• CQI Feedback: Each user feeds back his CQI if it is higher than ( on all shared channels. 
Otherwise, the user remains silent. 

• Compressive Sensing: BS finds the strong users using CS. 

• Least-squares estimation/refining: BS estimates or refines results obtained via CS using least- 
squares. 

• Optimum CQI Back off: BS backs off on the noisy CQI (SINR) based on the noise variance 
such that the throughput is maximized. 

3) User Selection: Select users and schedule them to beams. 

Remark 1: Once CQI has been determined for each beam direction, the base station can proceed to 
implement any of the various multiuser scheduling techniques. For example, the BS can go for random 
beamforming and might opt for a second-stage feedback to design the final precoding matrix [6]. The 
second-stage feedback is requested from the selected users only. Thus, the amount of second stage 
feedback is relatively much smaller compared to the amount of first-stage feedback. Alternatively, the BS 
could also implement the semi-orthogonal user selection (SUS) algorithm (zero beamforming) proposed 
in [5]. 

D. Feedback Protocol for the Digital Feedback Case 

The digital feedback is similar to analog feedback except that user feeds back "1" if his CQI for 
a particular beam is above a particular threshold. Otherwise, the user remains silent. To increase the 
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feedback granularity, we let the users compare his CQI to a set of thresholds, not just one. Thus, suppose 

that we want to set k thresholds (i < Q2 < , < Cfe such that the number of users whose CQI lie 

between the two consecutive thresholds [Q, Ci+i) is equal to s. Note that the last interval is 00). 
Following our discussion in subsection IIII-BI we can set the lowermost threshold as 



The feedback procedure is as follows: 

1) Threshold Determination: BS decides on thresholding levels Ci> C2? ••■> Cfc based on the sparsity 
level that can be recovered. For each threshold interval [Q, Q+i), repeat the User Feedback step. 

2) User Feedback: Repeat the following steps for each beam. 

• CQI Determination: Each user determines his best beam (corresponding to the highest CQI 
value). 

• CQI Feedback: Each user feeds back his CQI if it lies in threshold interval [Q, on all 
shared channels. Otherwise, the user remains silent. 

• Compressive Sensing: BS finds the strong users using Compressive Sensing. 

• Least-squares estimation/refining: BS estimates or refines results obtained via CS using least- 
squares. 

3) User Selection: For each beam, BS randomly selects one of strong users of the highest active 
threshold interval, where active threshold interval here means that there is at least one user sending 
feedback data in the interval. Here, CQI is the lower limit of the highest active threshold interval. 

IV. Throughput in the RBF Case 

In this section, we present the sum-rate throughput achieved by the proposed schemes. Although we 
focus on RBF, the proposed schemes can be applied to other beamforming methods (e.g. ZFBF). 

A. Throughput in the Analog Feedback Case 

The sum-rate throughput achieved in the RBF case with dedicated ideal feedback links is given by [4] 




Continuing in the same way, we get 





p 



K«E \ log 2 (l + max SINRj )TO ) 




(6) 
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Also, it is shown in [4] that ((6]) is equivalent to 

ft ~plog 2 (l + plog(n) -p{p- 2)loglog(n)) (7) 

As the SINRs fed back by the users are transmitted as is and the feedback links are noisy, so there is a 
need to back off the noisy received SINRs based on the noise variance as follows: 

SINR' = SINR + w (8) 

where the actual and noisy SINRs are denoted by SINR and SINR' respectively and w represents noise. 
Now, if we decide to back off the received SINRs by an amount A, then the back-off efficiency (rj) i.e. 
the probability that this backed off SINR is less than or equal to the actual SINR is given as follows: 

7] = P[SINR' - A < SINR] = ¥[w < A] = 1 - Q (J^j 

where Q represents the Q-function. Thus, the effective throughput (with back-off on noisy SINR) can be 
written as: 

ft e //~ (l-Q(^))plog 2 (/3- A) (9) 

where (3 = 1 + p\og{n) — p(p — 2) loglog(n). 

Differentiating ft e // w.r.t. A and setting it equal to 0, yields 

<£K^) exp (-£) log(/J - A) ^ (10) 

Simulation results confirm that the value of A that satisfies the above equation maximizes the effective 
throughput. 

B. Throughput in the Digital Feedback Case 

The sum-rate throughput achieved by p beams in the multiple thresholds (k in number) based digital 
feedback case for RBF is given below. 



K^pE 



log 2 (l + max Q 



where max Ci is the lower limit of the COI of the highest active threshold interval. 



July 10, 2009 



DRAFT 



11 



Alternatively, the same throughput can be derived analytically as follows. The throughput achieved for 

any transmit beam m is given as follows: 

k 

lZ m = ^log 2 (l + Ci)P(selected user in the threshold interval) P (threshold interval) 

i=i 

The probability of the threshold interval (denoted as Qi) is given by P(Qj) = [F(d + i) — F(Q)], where 
F(() is the cumulative distribution function (CDF) of CQI (SINR) defined as: F(Q = P[SINR < C] = 
1 — ^i+£) P - P i ; C > [4]. Capitalizing on the work of [15], we calculate the probability that selected user 
is in threshold interval Qi as follows: 

ra-l 



1 f n — 1\ 

P(selected user is in QA = V ( U ) V{P 2 

T^3 + 1 \ 3 J 



where 

V\ = P(j users other than the selected user are in Qi) = [F(£i+i) — F(Q)] 3 , and 
V 2 = P((n -3-1) users lies below the interval Q { ) = [F(Ci)] (n_j_1) 
Substituting these values of V\ and V2, and after some manipulations, one can show that 

P(selec«d user ,s m Qi ) = [F(0+l) _ F(0)| 

Thus, 

fc 

TZ rn = J2 log 2 (l + G)([^(Cm)] n - [^(0)] n ) 
1=1 

As, in our case there are p beams and all of them are identical, so the sum-rate throughput is given as 

p k 

TZ=J2 n m =pY J log 2 (l + Q)([F(Q +1 )} n - [F(C,)n 

m=l i=l 

V. Performance Evaluation 
We consider following metrics for the performance evaluation of the proposed feedback schemes. 

A. Feedback Resources Reduction 

There is a significant reduction in number of feedback channels required for carrying feedback infor- 
mation. The proposed schemes requires only 0(log(n)) feedback channels (shown in the Lemma given 
below) as opposed to n feedback channels required in the dedicated feedback case. 

Lemma 2: The number of multiple access feedback channels required for the proposed schemes is 
|(slog(n)), where c is a constant. 
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Proof: Specifically, let's assume that there are r channels shared between users over which feedback can 
take place. We can represent these channels using the system of equations (0. As already mentioned, 
(f2]) is similar to ones considered in [13]- [14], except that in our case the measurement matrix A, and 
the noise vector w are complex instead of real. So, we replace the complex-valued model in (O by its 
real-valued equivalent as shown below 

R(A) -9(A) 
9(A) 3fc(A) 

where 3i(A) & 9(A) represents real and imaginary part of A. After simplification the above equation 
reduces to 
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or, 



y = A v + w 



(11) 



The entries of A are i.i.d. M (0,1/2), and the entries of w are i.i.d. M (0,cr 2 /2). The above model 
dV-A| ) gives us the 2r x n real measurement matrix, and 2r x 1 real noise vector, so the sparsity pattern 
recovery techniques discussed in Section IIII-AI can be applied. Also, note that small values of s are 
sufficient (Section IIII-BI ). Therefore, we have 



2r « cs log(n), 



r = 2( sl °s( n ))- 



(12) 



Lemma 3: In the RBF case when n — > oo, the minimum number of multiple access feedback channels 
required is (log log log(re)) log(n). 

Proof: From Lemma 2, we have r = |(slog(n)) and for n — > oo, c = 2 [13]. For RBF systems with 
large number of users (n — > oo), the minimum value of s (the number of users who should feedback) 
required to achieve the sum-rate throughput is given by logloglog(n) [16]. Substituting these value of 
c and s in r = |slog(n), we conclude that number of multiple access (shared) feedback channels is 

r = (log log log(n)) log(n). 

B. Feedback Noise Reduction in the Analog Feedback Case 

The other important benefit of this scheme is the feedback noise reduction (which eventually results in 
better throughput) in the analog feedback case. This is because the feedback data of each user is carried 
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over all shared channels. This come in contrast to dedicated channel feedback case where feedback 
is carried over one channel only. We analyze the error covariance matrix (ECM), as it will allow us 
to identify the optimum amount of back off required on the noisy SINR which depends on the noise 
variance. Also, we analyze two measures of ECM — trace and determinant of error covariance matrix. 

1 ) Shared Feedback Channels: Error covariance matrix after the sparsity pattern is identified and LS 
is applied is given by [17] 

ECM = [Rv" 1 + AjRw^As]" 1 (13) 

where R v = Efvsvlj] = cr^I, and R w = E[w£Wg] = and ws refers to the entries of v and w 

corresponding to S. Substituting these values in ( [TBI , yields 



E[ECM] = E A - a 

(a) 



1 T 1 A* * 



(b) a , 



+ r 



(14) 



(15) 



where (a) follows because for fixed s and large r, E[A^As] — > rl [18], and (6) follows because for 
large r and high SNR (^f + r) —> r. For o\ = \, o%, = \j p where p is SNR. 
. Trace of ECM: Using (fl4l . 



E[tr(ECM)] = E|| ai || 2 ..|| aJa 



2\-l 



i=l 



sE|| a || 2 [(l + p||a| 



\2\-li 



where ||a|| 2 is a chi-squared variable with 2r degrees of freedom. Therefore, 

f-OO 1 

E[tr(ECM)l = s / (1 + px)- l —-x^ r -^e- x dx 

Jo r(r) 

where T(-) is the gamma function [19]. Using (3.383.10) of [19], the above evaluates to 

E[tr(ECM)] = s{l/p) r e 1/p T{\ - r, l/p) 

where T(-, •) is the incomplete gamma function [19]. 
• Determinant of ECM: Using (fl4l . 



E[det(ECM)] = Ea* a s 



det 



I+4rA^A 5 



(j- 
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We can write the above equation in terms of eigenvalues of A* S A$ as follows 



E[det(ECM)] = E V . A , JJ ( -L + 4) 




=rW(rt)" 

i0 °"w/ 



p(A) is given by ± E?=o 7i^b)T[ L r^)] 2 ^~ S e~ A [18], where If-'(A) = jj-e^A s_r- ^ r (e~ A A r ~' s+ ') 



is the associated Laguerre polynomial of order I. 
Alternatively, if we use the approximation in (031 ). trace and determinant of ECM are given by and 
(^■^p^j respectively 
2) Dedicated Feedback Channels: 



For a% = 1, a%, = 1/p where p is SNR. 

. Trace of ECM: Using dJ), tr(ECM) = n (1 + p) -1 . 

. Determinant of ECM: Using OH), det(ECM) = (1 + p)' n . 

3) Comparison between Shared & Dedicated Feedback Channel Cases: Trace of ECM and its ap- 
proximation is plotted in Fig. Q] Similar plot for determinant is omitted due to space limitation. Thus, 
from (031 ) & (fT6l) . we conclude that the back off on the SINR is 0(^=) in the shared feedback channels 
case as opposed to 0(<r w ) in the dedicated feedback channel case. 

Another point that needs to be noted is that the trace of ECM (which is commonly refereed as cost 
function and should be minimized [17]) in the shared feedback channel case is much smaller than that 
obtained in the dedicated feedback channel case. 

C. Feedback Load Reduction 

In addition to the feedback resources reduction, there is a reduction in the amount of feedback. In 
RBF scheme with dedicated feedback channel [4], n real values and n integer values (nlog 2 p bits) are 
fedback, as there are n users in the system. 

1) Analog Feedback Case: The proposed CS based analog feedback scheme requires only pr real 
values to be fedback. This is because there are r shared channels and the scheme is repeated for each 
beam. Note that the feedback load reduction is more dominant in systems with large number of users, 
as r ~ 0(log(n)) and p is small. 




(16) 
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2) Digital Feedback Case: The proposed CS based digital feedback scheme requires only pkr bits 
to be fedback. This is because there are r shared channels and the scheme is repeated for each beam 
& threshold. Note that the feedback load reduction is more dominant in systems with large number of 
users, as r ~ 0(log(n)) and p & k are small. 

D. Trade-off in the Digital Feedback Case 

Given a budget of bits that can be fedback, using intuition, it was shown in [20] that trade-off exists 
between the multi-user diversity and feedback accuracy. In our context, multi-user diversity is related to 
the number of shared channels r whereas feedback accuracy is related to the number of thresholds k, 
and so a similar trade-off may exist. The number of shared channels and thresholds must be chosen such 
that the throughput is maximized. This is explored using simulation in section IVIIi 

VI. Feedback Channel Training 

In the previous sections, we assumed that the channel A estimation is given to the system with the aid 
of a "genie" at no cost. In this section, we present how the feedback channel training can be accomplished 
and explore ways to reduce it. Here, we assume that $2} represents frequency feedback channels i.e., the 
entries of A, a,ij represents the gain of the j-th user in the i-th frequency band. 

A. Channel Matrix is Full 

The optimal number of symbols required for channel training is equal to the number of transmit 
antennas [21]. So we need p training symbols for the downlink channel and n training symbol for the 
uplink channel (as there are n users each having one transmit antenna). Training for each user in the 
uplink can be performed one by one, i.e., the first symbol of the coherence interval is reserved for user 1 
to perform training for all shared channels, and second symbol reserved for user 2, and so on. Continuing 
in this way, we need n symbol time to accomplish training for all users. Also, it is important to note 
that as there is little data to be sent for feedback purposes, so much of the uplink coherence time can be 
used for feedback training. Coherence time is typically of the order of few thousand symbols, so training 
would not be an issue for systems with moderate number of users. However, a method for reducing the 
amount of feedback channel training time is discussed in the next subsection. 

B. Channel Matrix is Block Diagonal 

In order to reduce the feedback training time, we divide the users into groups with each group being 
allowed to feedback only on a set of feedback channels, thereby reducing the full channel matrix to a 
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block diagonal one 



Ai 



A B d 



A fc 



Compressive sensing is applied in the same way as discussed in Section [TTlJ the only difference being 
that it is now applied on each block. Strong users in each block (or group) are found and the user 
corresponding to the maximum SINR among the strong users from all blocks is selected. As the users 
are i.i.d., so we divide the feedback resource equally among the k groups. Thus, training can now be 
performed for each block simultaneously. This approach reduces the feedback training time considerably, 
e.g. if we divide the total number of users into two groups, then the training will require n/2 symbol 
time as opposed to n symbol time required for the case when the channel matrix is full. 

The flip side of this approach is that compressive sensing is now applied on the group of users instead 
of all users as one block. Thus, for same sparsity level s (overall), with block diagonalization, the number 
of feedback channels required is given below 



Note that from the above equation it may first appear that the number of channels have reduced as the 
quantity inside the logarithm is reduced by a factor of k, however, it is the other way round. This is 
because now d has increased as the problem dimension (n) is reduced by a factor of k [13]. Thus, there 
is a trade-off between the reduction in the amount of feedback training and the number of feedback 
channels. Also, note that there is now an additional constraint requiring s/k to be an integer. 

C. Non-fading Channels 

When the channels are non-fading i.e., the channel gains are constant (or 1), then each strong users 
multiples its CQI with a unique binary chip sequence (consisting of ±1 each with probability 0.5) of 
length equal to the number of shared feedback channels r and send it over the multiple access shared 
channels. There are two ways of assigning chip sequences to the users: pre-programmed in users' device 
or sending it over the air. If it is send over the air, then the training time (used in the case of fading 
channels) can be used to send unique binary chip sequences to all n users. Thus, A is r x n Bernoulli 
matrix and so CS can be applied as Bernoulli matrices are shown to satisfy the RIP [10]. 



/a 



■BD 



Ia 1 -\ h /a* 
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VII. Numerical Results 

In this section, we present numerical results for CS -based feedback schemes by applying it in RBF 
context. We use p = 4 base station antennas, and n = 100 users. We set the threshold according to the 
sparsity level s, and use the maximum correlation technique (unless mentioned otherwise) for compressive 
sensing as this is much more computationally efficient than LASSO. Each point in the figures represents 
the sum-rate throughput achieved for shared number of channels determined by c & s according to (fl~2l) . 
We use SNR = 10 dB for both downlink and feedback link (unless stated otherwise) for calculating the 
sum-rate throughput. 

A. Analog Feedback Case 

In Fig. [U we present the sum-rate throughput with shared channel feedback in the analog feedback 
case. We use optimum back off on noisy SINRs in the analog feedback case. From this figure, we note 
that for small values of s the throughput is low. This is because the threshold works well for systems with 
large number of users but for systems with moderate number of users, we may have more or less number 
of users above the threshold than desired. So, if we set s low, then the probability that a beam has no 
strong user is relatively higher (resulting in a multiplexing loss) to the case when s is large. However, 
large values of s requires more feedback channels. Also, we see that the number of shared channels 
required to achieve the maximum possible throughput obtained in a noisy dedicated feedback scenario 
is 11 (corresponds to c/2 = 0.4 and s = 6). Also, it worth mentioning that the proposed scheme comes 
close to achieving the throughput obtained in a noiseless dedicated feedback scenario (dedicated feedback 
with ideal feedback links) due to feedback noise reduction. Note that 90% of throughput in noiseless 
dedicated feedback case is achieved by 19 shared channels (corresponds to c/2 = 0.8 and s = 5). 

In Fig. [3j we present results on block diagonalization method proposed for reducing the feedback 
training time (section IVI-Bt . Here, we divide 100 users into two groups of 50 users and compressive 
sensing is applied on each group. It is clear from the figure that this method requires few more feedback 
channels. Also in this figure, we present result based on LASSO which shows that LASSO method 
performs marginally better than maximum correlation method. 

In Fig. @1 we present the sum-rate throughput achieved by two-stage RBF in the analog feedback case 
when the feedback channels are noiseless. In the second-stage of the two-stage RBF, additional feedback 
information (beam gain information (BGI) [6]) is requested from the selected users only and Iterative 
Beam Power Control (IBPC) algorithm proposed in [6] is used for the re-distribution of the total power 
among the active beams (beams for which there are strong users) in an optimized manner. From the figure, 
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we see that there is hardly any gain for two stage RBF with dedicated channel feedback, however, it is 
evident that for compressive sensing based opportunistic feedback protocol, two-stage RBF is effective 
even for moderate to large number of users. This is because if no user is strong for some beams, the 
system still suffers from the multiplexing loss but the power of those beam are distributed among the 
active beams in an optimized way Also, note that there is no back-off required here as the feedback links 
are noiseless. With two-stage RBF, the number of shared feedback channels required is 15 (corresponds 
to c/2 = 0.8 and s = 4). 

B. Digital Feedback Case 

For all digital feedback cases, we chose s = 1 (the minimum possible value) and set multiple thresholds 
as discussed in section IIII-Di This is because for the proposed scheme, s = 1 will allow us to set the 
highest possible uppermost threshold thereby ensuring a higher throughput. 

In Fig. [5J we present the sum-rate throughput achieved with shared channel feedback in the digital 
feedback case. Is is evident form the figure that the proposed scheme in a noisy scenario achieves the 
throughput obtained in a noiseless dedicated feedback scenario (dedicated feedback with ideal feedback 
links). Also, we see that the throughput increases with the increase in the number of shared channels & 
thresholds. Taking the pessimistic view, we need only 10 feedback channels (corresponds to c/2 = 2 and 
s = 1). However, it is important to note that beyond a certain number of shared channels or thresholds, 
the throughput either becomes stagnant or increases marginally. 

In Fig. [6j we consider fixed budgets of p x kr bits that can be fedback. From the figure, we note that 
such a trade-off exists and for a given fixed budget there is an optimum number of thresholds and shared 
feedback channels that maximizes the throughput. 

VIII. Conclusions 

In this paper, a generic feedback channel model and compressive sensing based opportunistic feedback 
schemes are proposed. The proposed generic feedback channel model is shown to encompass all existing 
feedback channel models proposed in the literature. We have shown that the proposed analog & digital 
opportunistic feedback schemes achieves the same sum-rate throughput as that achieved by dedicated 
feedback schemes, but with feedback channels growing only logarithmically with number of users. Also, 
we derived an expression for the sum-rate throughput in the digital feedback case with multiple thresholds. 

In the analog feedback case (noisy scenario), it has also been shown that due to feedback noise 
reduction, the proposed scheme comes close to achieving the throughput obtained in the case of noiseless 
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dedicated feedback. In the digital feedback case, it has also been shown that beyond a certain number 
of shared channels or thresholds, the throughput either becomes stagnant or increases marginally. Also, 
given a budget on the amount of bits that can be fedback, we have shown that there exist a trade-off 
between the number of shared channels and thresholds and therefore they must be chosen such that the 
throughput is maximized. 

Although the results presented here only show the performance of the the proposed schemes in the 
RBF context, the schemes can easily work with other beamforming methods. 

IX. Appendix 
Compressive Sensing 

Here, we give the reader a brief introduction about compressive sensing. Let v g K" be an unknown 
vector, with at most s non-zero entries (s < n/2) and let S denote its support set with |5| = s << n. 
Suppose that we make a set {yi, . . . , y r } of r independent and identically distributed (i.i.d.) observations 
of the unknown vector v, each of the form 

yi = af v + Wi (17) 

where wi ~ Af (0, N ) is observation noise, and a» ~ Af (0, I nX n) is a measurement vector. In the matrix 
form, it can be compactly written as 

y = Av + w (18) 

Reconstruction will not be possible if the measurement process damages the information in v, which 
often happens in practice. A necessary and sufficient condition for the system of equations to be well- 
conditioned (thus having a stable inverse) is the restricted isometry property (RIP) [10]- [12]. 
Definition 1: A r x n matrix A has the s-RIP with appropriately chosen constant < e s < 1 if 

||Av||» 

1 - e s < <l + e s (19) 

IMIa 

holds for all s-sparse vectors v. 

The s-RIP property ensures that the matrix A preserves the lengths of these particular s-sparse vectors. 
Practical recovery algorithms require that A satisfies a more conservative RIP (3s-RIP in general [12]). 
If A is a random matrix consisting of Gaussian random variables, the RIP property is satisfied with 
overwhelming probability [11]. In this case, the number of measurements that are necessary to recover 
v efficiently in a noiseless scenario with high probability is on the order of r ~ slog(n/s). 
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Figure 1 . Trace of Error Covariance Matrix for shared feedback 
channel (SFC) and dedicated feedback channel (DFC), n = 100 
and uplink SNR = 10 dB for different values of s. 



Figure 2. Analog Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink 
& feedback link) for different values of s. 
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Figure 3. Analog Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both 
downlink & feedback link) for different methods viz. LASSO, 
maximum correlation, and maximum correlation with block 
diagonalization. 



Figure 4. Analog Shared Channel Feedback: Throughput versus 
c/2 for Single stage and Two stage RBF, p = 4, n = 100 and 
SNR = 10 dB & oo for downlink and feedback link respectively. 
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Figure 5. Digital Shared Channel Feedback: Throughput versus Figure 6. Digital Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink k for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink & 
& feedback link) for different values of k. feedback link) for different budgets of bits that are to be fedback. 
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